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Abstract. Dirac's constraint analysis and the symplectic structure of geodesic 
Q^ ' equations are obtained for the general cylindrically symmetric stationary spacetime. 

For this metric, using the obtained first order Lagrangian, the geodesic equations 
of motion are integrated, and found some solutions for Lewis, Levi-Civita, and Van 
^ ■ Stockum spacetimes. 
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1. Introduction 

Dirac first worked out tlie tlieory of quantizing constrained systems in general|lj, and 
^ ■ general relativity in particular and his pioneering work continues to serve as the 
^ • foundation of current efforts to canonically quantize gravity. Besides the hamiltonian 
formalism utilizing the Schrodinger representation, there is an alternative Hamiltonian 
approach which is based on Dirac's analysislHj of constrained systems. Many more 
extensive studies of the subject can be found in the literature; see, for example, Refs. 4 
and 5. In this context, we discuss Dirac analysis of geodesic equations for the general 
cylindrically symmetric stationary spacetimes, and later integrate the obtained first 
order equations of motion for these spacetimes. 

Because of both the mathematical simplicity and the physical relevance to our 
realistic worldjHl Ej, space-times with cylindrical symmetry have been extensively 
studied, and their relativistic applications have been further discussed recently |H]~[I21- 
The general form of this metric in vacuum case was given by Lewis^S]- Lewis stationary 
vacuum metric is usually presented with four parameters ^3] which admits a specific 
physical interpretation when matched to a particular source. These four parameters 
which are related to topological defects [HI El not entering into the expression of the 
physical components of curvature tensor may be real (Weyl class) or complex (Lewis 
class). In recent years, the physical meaning of these parameters have been discussed for 
both classes jni HDj- The corresponding static limit of the the Lewis class was obtained 
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by Levi-Civita (LC)Pni- Even in the simplest case of the LC solution, its physical 
interpretation is not completely understood, yet. In general, it contains two independent 
parameters, in which there is only one mass parameter. One of them is associated with 
the topological defects while the second parameter is connected with the mass per unit 
length. Another special case of Lewis metric is the van Stockum solution ^2] which 
represents the gravitational field produced by a rigidly rotating dust cylinder with a 
finite thickness. The matching of this space-time to the vacuum Lewis space-time was 
also completed in Ref.l2, and studied in detail by Bonnor 18J. 

The paper is organized as follows. In the next section we present the stationary 
cylindrical spacetime in general and give the familiar spacetimes at the exterior of the 
boundary of the source. In Sec. El we shall present the symplectic structure of geodesic 
equations for the general cylindrically symmetric stationary metric. For the first order 
form of the Lagrangian that yields the geodesic equations for this metric, we shall apply 
Dirac's theory of constraints [3] to the degenerate the Lagrangian. We find that the 
constraints are second class as in the case of all integrable systems. The constraint 
analysis yields the Dirac brackets, or the Hamiltonian operators in the language of 
integrable systems. The symplectic 2-form is obtained by the Poison bracket of Dirac's 
constraints which is also the inverse of the Hamiltonian operator. In Sec. |31 using 
the first order Lagrangian for the general cylindrically symmetric metric, the geodesic 
equations are integrated, and the obtained results are worked in cases of the Lewis 
spacetime for the Weyl class, the exterior van Stockum and LC spacetimes. 

2. Spacetime 

The general line element for a cylindrically symmetric stationary spacetime is given by 

ds^ = -fdf + 2kdtd(j) + e'^((ir2 + dz^) + (1) 

where f,k,n and i are functions only of r, and = {t,r, z,(f),t),i = 0, 1,2,3 are the 
usual cylindrical coordinates with 

-oo<t,2<oo, r>0, 0<(/)<27r (2) 

and the hypersurfaces (p = 0,2n being identified. Einstein's field equations for vacuum 
are 

R^J = 0. (3) 

The general solution of (jH)) for is the stationary Lewis metric which can be 
written as 



/ = ar-"+i - ^r"+\ (4) 



k = - Af, (5) 
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with 



A = + h. (8) 

naj 

The constants r?,, a, h and c can be either real or complex, the corresponding solutions 
belong to the Weyl or Lewis classes, respectively. For the Weyl class, the above 
parameters have the following physical interpretations. The parameter n is associated 
with the Newtonian mass per unit length of an uniform line mass a when it produces 
the low density regime. The parameter a is related to the constant arbitrary potential 
that exist in the corresponding Newtonian solution, while the parameters h and c are 
responsible for the non-staticity of the spacetime, since when we take 6 = and c = 
the Weyl class reduces to the static LC metric. The parameter h is related, in locally 
flat limit, with the angular momentum of a spinning string. The parameter c measures 
the vorticity of the source when it represented by a stationary completely anisotropic 
fluid. For further details see Ref. 9. For the Lewis class, the physical and geometrical 
meaning of the four parameters of the Lewis metric was given in Ref. 10. 

For the line element of the LC static vacuum spacetime in the Weyl form, the 
functions /, k, /i, and £ have the following expressions [HI 13 HI] 

f = r^- k = 0, i = C-^r^-"'' , e^ = r*-{^-^) (g) 

where a and C are two arbitrary constants and both of them are fixed by the internal 
composition of the physical source. The constant C refers to the angular defect [H], and 
cannot be removed by scale transformation. This constant is related, in the locally fiat 
limit, with the parameter a in the Lewis solution for the Weyl class, given by 

a = C^. (10) 

The physical importance of the other parameter a is mostly understood in accordance 
with the Newtonian analogy of the LC solution, i.e. the parameter a represents the 
mass per unit length jHl EHl- The parameter a is connected to the parameter n in the 
Lewis spacetime for the Weyl class as 

n = l-4a. (11) 

The functions /, k, /i, and £ for the van Stockum solution^Jj are given by 

/ = 1, fc = ar^ £ = r2(l-aV), ^ = -aV^ (12) 

with a being an arbitrary positive constant. The energy density and the four velocity 
of the dust are 

ZnU 

where G is the gravitational constant. The angular velocity to the fluid with respect to 
a locally nonrotating frame is = a{l — a^r^)"^. Since near the axis, a; ^ a, it can be 
interpreted that a is the angular velocity of the fluid on the axis|18j. 

The van Stockum exterior solution (fT^ . which is a particular case of the Lewis 
metric, contains the globally Minkowski spacetime as a special case. Therefore the 
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van Stockum solution (fT^ must be a particular case of the Weyl classj^]. Since the 
van Stockum spacetime cannot be reduced to the globally static LC metric ©, it is a 
particular case of the Weyl class with b ^ and c 7^ 0, since for 6 = and c = the 
Weyl class can be globally reduced to the static LC metric 



3. Dirac Analysis for Geodesic Equations 

The equations governing the geodesies can be derived from the lagrangian 

where r is an afiine parameter along the geodesies. From the external problem it emerges 
the Euler-Lagrange equations 
d dC 
dr \dxy dx^ 
and from them follow the geodesies given by 









(14) 



/Vi^ _|_ 



(15) 



where the overdot denotes differentiation with respect to r. For spacetime the 
Lagrangian (fT^ is 

O \ 1 o 



(16) 



This Lagrangian is second order and therefore not suitable to a discussion of symplectic 
structure. For purposes of Hamiltonian analysis we need to start with first order 
Lagrangian and it can be verified that 



k{P + Q) - D{P - Q) 



t 



1 



iP + 



- -e^f^iM^ + N^) + Mr + Nz 
gives rise to the equations of motion 

P = ^+{k + D)i, Q = e^+{k-D)i, 



k{P + Q) - D{P - Q) 

i 



0, 



[P-Q) 

AD 



P + Q = 0, 
iV = 0, 

which together result equations ()A.1IIA.3]) in Appendix A[ where we have defined 

D^ = k^ + if 



(17) 

(18) 
(19) 

(20) 

(21) 

(22) 
(23) 
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and 



F^{P + Q)[-] -{P-Q){- 



D 



and the prime represents derivative with respect to r. Dirac quantization is canonically 
quantize the original phase space which is usually even dimensional symplectic manifold, 
and then imposed the gauge constraints as operator conditions on the physical quantum 
states. In this first order formulation we have introduced x"' = {P.Q, M, N), a = 
4, 5, 6, 7, as new variables which is double the number required. So we consider a 
symplectic manifold spanned by variables = {x\x"'),A = 0,...,7, where a;*'s and 
x^'s are, respectively, spacetime and configuration space variables. Then, first order 
field equations become 

= X(X^), 

with the vector field defining the flow 



(24) 



X 



1 

2b' 



+ 



+ 



P + Q-^{P-Q) 



6 6 
M— + N— 
or oz 



2i 



6 Mie-^" 



-e-^f^'{M' + N') + —PQ + 



F 
2D 

{P-Q) 

AD 
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SM 



(25) 



det 



for the geodesic equation (fT3j) . 

The Lagrangian (jl7j) is degenerate because its Hessian 

. . =0 (26) 

dx^ dx^ 

vanishes identically. Hence it is a system subject to constraints and the passage to its 
Hamiltonian structure requires the use of Dirac's theory of constraints^. We introduce 
the canonical momenta of the test particle defined by 

dC 



H. 



dX^ 



(27) 



which cannot be inverted due to equation The definition of the momenta therefore 
gives rise to the constraints 



$0 

$1 

$2 
$3 

$5 
$6 
$7 



Hi 
H, 



H,-iV, 

^^-\{P + Q). 

Hp, 
Ha/, 



(28) 
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which must vanish weakly, i.e on shell. In order to determine the class of these 
constraints we need to obtain the Poisson bracket of the constraints 



CAB(r,T) 



$A(r),$B(r^ 



using the canonical Poisson brackets 



(29) 
(30) 



between the dynamical variables and their conjugate momenta. The result 
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5(r-f)(31) 



shows that the constraints (j28jl are second class as in the case of all integrable 
systems P3J. 

In order to obtain the Hamiltonian for the degenerate Lagrangian (jl7|) we first 
construct the free Hamiltonian obtained by Legendre transformation 



(32) 



and the total Hamiltonian density of Dirac is given by 

Ht = Ho + X^^a (33) 

where are Lagrange multipliers. Since we have second class constraints the Lagrange 
multipliers will be determined from the solution of 

Ht,<^a 



(34) 

which ensure that the constraints hold for all values of r. Since the constraints are linear 
in the momenta the Lagrange multipliers are given by 











= Me"^, 




A^ 






A^ 




A^ 


2D ' 




A^ 


= -A^ 




A^ 


i' 1 


'y(M' + N') + ^{P - Q), 


A^ 


= 





(35) 



7 



which follows directly from the flow (|25|) . 

The Dirac bracket is a modification of the Poisson bracket designed to vanish on 
the surface defined by the constraints. For two smooth functionals A, B of the canonical 
variables we have 

{A, B}d = {A, B} - {A, <l>^} J^^{<I>B, B} (36) 
where J is obtained by inverting the matrix of the Poisson bracket of the constraints C 

Cab^t, f ) J^^(f , f ) (B = b{T - f ) (37) 

The inverse of the Poisson bracket of the constraints is known as the Hamiltonian 
operator in the literature of integrable systems [22j. From Eq. (j37j) we obtain 
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and Eqs.()18 |) - (f^ can be written in Hamiltonian form 



J 



AB 



5X^ 



\ 


1 





(39) 



5(r-r)(38) 



where integration over dotted variables is implied. 
The symplectic 2-form is given by [221 
1 



B 



(40) 



and using Eq.(|^ we find that 



-SPA 
2 



+ -SQA 



+ 5M A5r + 5N A5z + 



JJ 



5P A 



(42) 

of the symplectic 



+ ^-^if^^Q A 5t + i [(P + Q) (j)' -{P-Q) i^)] Sr A St (41) 
Finally, Hamilton's equations can be written in the form 

where ix denotes contraction with respect to the vector field 
2-form dH}. 

We shall now turn to the Witten-Zuckerman|2()| E] formulation of symplectic 2- 
form vector density uj which is closed 

6u = (43) 

and conserved 

u = 0. (44) 
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Starting with the Lagrangian (jl7|) we find that the Witten-Zuckerman symplectic 2-form 
is the same expression for the symplectic 2-form ()41|) obtained from Dirac's theory of 
constraints, i.e. 



4. Integration 



The geodesic equations (fT3j) are four equations for the four unknowns i, r, i and (j) (see 
[Appendix A] ). Firstly, from the Eqs. (|T8|) and (|T9|l . we obtain 



i = — (P-Q), 



r 

z 



1 



2i 



P + Q)--{P-Q) 



Later, integrating Eqs. PT|) -(| ^ we get 

N = K, 

P + Q =2L, 

j{P-Q)-j{P + Q)=2E, 



(45) 

(46) 
(47) 

(48) 

(49) 
(50) 

(51) 



where K, L, and E are integration constants. Using these results in Eq. ()2()j) and 
rearranging, yields the following Riccati type differential equation 

M = p{t)M^ + g(r) (52) 

where p and q are given by 



1 



2if 



{Lk + Ei) 



2D^ 



where F = 2 



{hf-Ek^ ~E^D 
^ il)' " ^^^g^^'' ("f )' • Then, the solution of Eq. (jSSI) is obtained as 

1/2 



(53) 
(54) 



M 



U 



where the following integral equation must still be satisfied for U, 



dU 



+ Ci 



p 



1/2 



dT 



U^ + PU+l J \q 

where Ci is an integration constant, and /? is given by 

2p\p J \ q 
Using (jlUj) in (jHSj), we find 
r = e'^'M. 



(55) 

(56) 

(57) 
(58) 
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Thus, we have obtained the generic expression for the radial speed of the test particle. 
A detailed examination of the solutions of the Eq. (j^^ is beyond the scope of this 
paper;therefore, we shall only consider some special orbits here. Let us first assume 
that p = —q = —1. In this case, we find the following solutions: 



(59) 
(60) 



e ^ = 2r + m, 

{1 ^^2g±2(r-ro) _ ^ /or M = ±1, 

I [n'^ cosh^(r — tq) — m] , for M = tanh(T — tq) 
I [n^ sinh^(r — Tq) — m] , for M = coth(r — Tq) 

where m, n, and tq are integration constants. Second, we assume that p = q = br"-, 
where b and n are nonzero constants. In this case, it follows from Eqs. ()52j) and ()53j) 
that for n 7^ — 1 



a cos 



To 



r r 



+ 1 

ab sign{b) cos(2To)r^+"' 



M = tan 



n + 1 



(61) 



'1 + n)(2 + n)sign{l + n 
1 



X 1F2 

osin(2ro)T 



.1 

^2 ^ 2(1 + n) 



3 3 
^'^2' 2 



2 1 + n 



^2 ^2(l+n) 



1-^2 



1 1 1 

^2(l + n)^'^2' ^2(l + n) 



1 + n)2 

f,2 ^2{l+n) 



:i + n) 



and for n = — 1 



a cos^ [feln(ror)] , M = tan[bln(ror)]. 



r(T) 



a r 



■cos[261n(ror)] [tan[261n(ror)] - 26] , 



(62) 

(63) 
(64) 



2(1 + 46)2 

where o and tq are constants of integration, and 1F2 is the generalized hypergeometric 
function. 

Considering the Eq. (j27j) . the momenta of the test particle for metric (pj using the 
Lagrangian (fTTj) is given by 

= M, 



n^^^(p + g) = L, 



(65) 



21 



[k{P + Q) - D{P - Q)] = -E 



Hence E can be interpreted as the total energy of the particle, and will be always 
taken nonnegative. K can be interpreted as its momentum along z and L its angular 
momentum. In terms of these conserved quantities, t, f, z and (j) become 

-{Ei + Lk) 



1/2 



(66) 
(67) 
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z = Ke-f", (68) 
• Ek-Lf 

0=— — • (69) 

As is well known, the Eq. has a first integral that is equivalent to g^^pX^x'^ = — e, 
where e = 0, 1, or —1 if the geodesies are respectively null, timelike or spacelike. This 
implies that 

-e = -fi^ + 2ki(j) + e^{r^ + z"^) +e<j)'^ (70) 

so that Cl in ffTBj) is equal to | along the path. For the Lagrangian (fTTj). this equation 
becomes 

jPQ + e-''{M^ + N^) = -e. (71) 

Now, substituting (jlKjl . (jTfj) and (jlKjl into (jHH), we have an expression for the radial 
speed of the test particle : 

f2 = e"'^ [Wo{r) - W{r)] , (72) 

where Wo(r) and W{r) are defined as 

Wo{r) = E'^ + 2EL^ - e, W{r) = L^i^ + K'e->^. (73) 
From the Eqs. (j67j) and (j72j) . we find 



U = ie'^/' [Woir) - W{r)] (74) 

which enables us to relate Eqs. ^7\i and (ff^ . On the other hand, setting Wo{r) = 
G{r) [Vo - Vi{r)] and W = G{r)V2ir), then Eq.(jZ2I) becomes 

f2 = e-f'Gir) [Vo - V{r)] , (75) 

where V = V1 + V2. Then, one can obtain the functions G{r) and V{r), and the constant 
Vo for spacetimes (HJ with dH)-©, © and ((121) given in Sec. El 

In the Lewis spacetime for the Weyl class, which means that the parameters n, a, b 
and c appearing in (jlll-Q are real, we obtain Vq, V{r) and G{r) as 

G{r) = r"-\ (76) 

Vo =^[{Eb + L)c-En]\ (77) 

V{r) = er^-" + a{Eb + Lfr'^'' + kV^^-^^^^^+'^Z^. (78) 

In the general case, i.e. 6 7^ and c o, we see that the parameter c only effects 
Vo by modifying the energy of the test particle, and leaving otherwise the geodesies 
indistinguishable from the static LC spacetime. In two different subcases, when 
6 = 0, c 7^ and 6 7^ 0, c = 0, the physical interpretations are given by Herrera and 
Santos [IS]. Furthermore, for the Lewis solution (jH)-©, the differential equation (f73|l 
can be written in the integral form as 

n-l 

/r ^ dr / X / N 

= ±(r-ro), (79) 
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where F = Eb + L, and tq is a constant of integration. The ± signs in Eq. (|79|) 
correspond to outgoing and ingoing geodesies, respectively. 

In the case of the LC solution Q, for Vq, V{r) and G{r), we have 

G{r) = r-^", (80) 
Vo = E\ (81) 
V[r) = e/'^ + c2lV2(4'^-i) + kV^'^^i-'^). (82) 

Then, it is seen that the parameter C, or a due to (fTUj). does appear in (f77j) . but does 
not in (jHH)- This is a different result obtained from the Lewis spacetime for the Weyl 
class taking 6 = and c = 0, in which the spacetime reduces to the static LC spacetime. 
For the van Stockum solution (fT^ . we find that Vq) V{r) and G{r) are as follows 

G(r) = 1, (83) 
^0 =E'^ + 2aEL - e, (84) 
y (r) = a^E'r^ + L^r"' + K^e'^"'" . (85) 

In this case, V{r) is non-negative and, in order to have Eq. (ffH|l meaningful for real r, 
we must have Vq > 0, which is equivalent to 

E> {a^L^ + ef/'^ -aL. (86) 

Finally, we note that Eqs. (jHEI), (lEHI), and (ffH|l describe the motion of test 

particles in the background of the general cylindrically symmetric stationary spacetimes. 
Using these equations, it can easily be obtained the orbits for the particles, and the 
effective radial potential. 

5. Conclusion 

In this paper, we have presented the Dirac analysis of geodesic equations of the general 
cylindrically symmetric stationary spacetimes in explicit form. Using Dirac's theory 
of constraints and the covariant Witten-Zuckerman approach we have obtained the 
Hamiltonian operators. The results for the symplectic 2-form coincide in both of these 
theories. We note that the original Lagrangian (fT^ which is second-order is non- 
degenerate, and gives second-order geodesic equations which are given in Appendix A 



However, the first-order Lagrangian ()17j) is degenerate and produces first-order geodesic 
equations. Therefore, if we consider this first-order Lagrangian given by (fT7j) . then 
we can easily find the solution of the obtained first-order geodesic equations of motion 
for the considered spacetimes. In the previous section. Sec. El using this degenerate 
Lagrangian approach, we have integrated the geodesic equations of motion for the 
general cylindrically symmetric stationary spacetimes, and found some solutions for 
Lewis, LC, and Van Stockum spacetimes. Also, for the radial speed of the test particle, 
we have found a generic expression given in (jKH|) which is depend on the solution of Eq. 
fl^ . In some spacial cases, we have solved the Eq. and found some solutions for 
the radial speed of the test particle. 
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Appendix A. Second-Order Geodesic Equations of Motion 

For the spacetime (0), it follows from the Lagrangian (fTT)|) that the geodesic equations 
of motion (fT3j) are 

Dt+^^ tr + f0 = O, (A.l) 

2r + e-^(/'t^ -2A;'t0-/>2) = 0, (A.2) 

•■ fk' - kf ■ fi' + kk' . , , , 

z + fi'rz = 0, D(j)+- — -jj-^f^ + —^ r(/) = 0. (A.3) 
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